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Abstract 

We present measurements of the temperature-dependent frequency shift of five niobium super- 
conducting coplanar waveguide microresonators with center strip widths ranging from 3 /jm to 
50 /jm, taken at temperatures in the range 100-800 mK, far below the 9.2 K transition temper- 
ature of niobium. These data agree well with the two-level system (TLS) theory. Fits to this 
theory provide information on the number of TLS that interact with each resonator geometry. The 
geometrical scaling indicates a surface distribution of TLS, and the data are consistent with a TLS 
surface layer thickness of order a few nm, as might be expected for a native oxide layer. 
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Superconducting microresonators have attracted substantial interest for low temper- 
ature detector applications due to the possibility of large-scale microwave frequency 



multiplexing 
puting experiments 10 
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9j. Such resonators are also being used in quantum com- 
and for sensing nanomechanical motion fl3|. We previously 
reported that excess frequency noise is universally observed in these resonators regardless 
of the type of superconductor or substrate being used, and suggested that two-level systems 
(TLS) in dielectric — may be responsible for this ooiseQ. TLS effect, are also 
observed in superconducting qubits[ll|. The TLS hypothesis is strongly supported by the 



observed temperature dependence of the noise, and also by the observation of temperature- 
dependent resonance frequency shifts that agree closely with TLS theory [15)]. To make 
further progress, it is essential to constrain the location of the TLS - to determine whether 
they exist in the bulk substrate or in surface layers, perhaps oxides on the exposed metal or 
substrate surfaces, or in the interface layers between the metal films and the substrate. In 
this paper, we provide direct experimental evidence for a surface distribution of TLS. 

The unusual physical properties of amorphous materials at low temperatures are con- 
sistent with a collection of TLS with a wide range of excitation energies and relaxation 



rates 
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171 ] . The TLS have electric dipole moments that couple to the electric field E of 



our resonators. For microwave frequencies, and at temperatures T between 100 mK and 1 K, 
the resonant interaction dominates over relaxation, which leads to a temperature-dependent 
variation of the dielectric constant given by [3] 
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where w is the frequency, \I/ is the complex digamma function, and 5 = nPd 2 /3e represents 
the TLS-induced dielectric loss tangent at T = for weak non-saturating fields. Here P and 
d are the two-level density of states and dipole moment, as introduced by Phillipsfisl]. 

Eq. (CQ) has been used extensively to derive values of Pd 2 in amorphous materials. If 
TLS are present in superconducting microresonators, their contribution to the dielectric 
constant described by Eq. <^ could be observable as a temperature-dependent shift in the 
resonance frequency. Indeed, it has recently been suggested that the small anomalous low- 
temperature frequency shifts often observed in superconducting microresonators may be 

and in fact excellent fits to the TLS theory can be obtained [ij]]. 



due to TLS effects 
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Assuming that the TLS are uniformly distributed in a volume Vh of host material (e.g. a 
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metal oxide or the bulk substrate) which has a dielectric constant of e^, it can be shown 
that the fractional resonance frequency shift is given by 

where the filling factor F is given by 

p _ Iv h ^Eiffdr _ K 
J v eE{r) 2 df w e ' 

The factor F accounts for the fact that the TLS host material volume Vh may only partially 
fill the resonator volume V, giving a reduced effect on the variation of resonance frequency. 
According to Eq. (J3]), F is the ratio of the electric energy w e h stored in the TLS-loaded 
volume to the total electric energy w e stored in the entire resonator. 

The key idea of the experiment described in this paper is to measure Af T / f r of coplanar 
waveguide (CPW) resonators with different geometries in order to obtain values of F5 for 
each geometry. The frequency-multiplexed resonators are all fabricated simultaneously and 
are integrated onto a single chip, and are measured in a single cooldown. We can therefore 
safely assume that a single value of the loss tangent 5 applies for all resonator geometries. 
This allows the variation of the filling factor F with geometry to be determined, providing 
information on the geometrical distribution of the TLS. If TLS are in the bulk substrate 
with dielectric constant e r , Eq. (j3J) applied to the CPW field distribution would yield a 
filling factor F ~ e r /(e r + 1) that is independent of the resonator's center strip width s r . 
If instead the TLS are in a surface layer, F should be dependent on the CPW geometry, 
scaling roughly as l/s r . 

We used a device with a 120 nm-thick Nb film deposited on a crystalline sapphire sub- 
strate, patterned into five CPW quarter-wavelength resonators with different geometries. 
Because Nb has a critical temperature T c = 9.2 K, the effect of superconductivity on the 
temperature dependence of the resonance frequency is negligible for T < 1 K. As shown in 
Fig. 1, each resonator is capacitively coupled to a common feedline using a CPW coupler 
of length l c = 200 /zm and with a common center-strip width of s c = 3 /mi. The coupler is 
then widened into the resonator body, with a center-strip width of s r = 3 /im, 5 //m, 10 /im, 
20 /jiii or 50 /tm, and a length of l r ~ 5 mm. The ratio between center strip width s and 
the gap g in both the coupler and the resonator body is fixed to 3:2, to maintain a constant 
impedance of Z « 50 Q. The resonance frequencies are f r ~ 6 GHz, and the coupler is 
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designed to have a coupling quality factor Q c ~ 50, 000. The device is cooled in a dilution 
refrigerator, and its microwave output is amplified using a cryogenic high electron mobil- 
ity transistor (HEMT) amplifier on the 4 K stage. The complex transmission S21 through 
the device and HEMT is measured using a vector network analyzer locked to a rubidium 
frequency standard, and resonance frequencies are obtained by fitting these dataj^, 15]. 

Fig. 2 shows the measured frequency shifts Af r /f r for the five resonators as a function 
of temperature over the range 100 mK to 800 mK. Although all of the resonators display 
a common shape for the variation of frequency with temperature, the magnitude of the 
effect varies strongly with geometry. As shown by the dashed lines in Fig. 2, fits to the 
TLS model (Eq. [2]) generally agree quite well with the data. The nonmonotonic variation of 
the dielectric constant with temperature is familiar from the TLS literature: f r increases (e 
decreases) when T > hw/2k; a minimum in f r (a maximum in e) occurs around T = hu)/2k; 
at lower temperatures (T < 100 mK), we would expect to see a decrease in f r (increase in 
e) as indicated by the extrapolation of the fit. The largest deviations from the TLS model 
(about 4%) occur at the lowest temperatures, and are likely due to TLS saturation effects. 
Indeed, power-dependent frequency shifts of this size are expected theoretically and have 
a ,o bee, p r ev 10U , y observed e„ ntally[1 i * this paper, we will ignore these small 
effects and focus on the geometrical dependence. 

With the exception of the 3 /im resonator, the measured values of F5 from the fits have 
to be corrected for the coupler because the coupler's center strip width s c = 3 /im differs 
from that of the resonator, s c 7^ s r . In the limit l c << l r , it can be shown that the corrected 
filling factor is given by F* = (F — £F 3/im ) /(l — t), where t = 2l c / {l c + l r ). The values of F*5 
are listed in Table 1, as well as the ratios relative to the value for 3 /im resonator. 

A portion of the CPW inductance per unit length is contributed by the kinetic inductance 
of the superconductor. The kinetic inductance fraction, a, depends on CPW geometry and 
may be determined by measuring resonance frequency shifts at higher temperatures, closer 
to T c [2jJ. We therefore measured the resonance frequencies at 4.2 K (0.46 T c ), allowing the 
shift A/ r (4.2 K) = / r (4.2 K) — / r (100 mK) as well as the kinetic inductance fraction to be 
calculated for each geometry, as shown in Table 1. 

Fig. 3 shows the results for the geometrical scaling of the corrected filling factor F* and 
the kinetic inductance fraction a, plotted as ratios relative to their respective values for the 
resonator with a 3 /im wide center strip. The observed strong variation of F* with geome- 
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try immediately rules out a volume TLS distribution, and favors a surface distribution. We 
investigate this in more detail by comparing the data to two theoretically calculated geomet- 
rical factors g m and g g , which have units of inverse length and are calculated from contour 
integrals in a cross-sectional plane given by g m = j mctal E 2 dl/V 2 and g g = j^ p E 2 dl/V 2 , 
where V is the CPW voltage. The first integral is actually a sum of three contour integrals, 
taken over the surfaces of the three metal conductors, the center strip and the two ground 
planes. The second contour integral is taken over the two "gaps" , the surface of the exposed 
substrate in between the conductors. These contours are illustrated in the inset of Fig. 1. 
The integrals are evaluated numerically using the electric field derived from a numerical 
conformal mapping solution to the Laplace equation. 

According to Eq. 0, F* should have the same scaling as g m if the TLS are distributed 
on the metal surface (or at the metal-substrate interface), or as g g if the TLS are located 
on the surface of the exposed substrate. The kinetic inductance of the CPW may also be 
calculated using a contour integral similar to that of g m , except that the integrand is replaced 
by H 2 21]. Because the magnetic field H is proportional to E for a quasi-TEM mode, we 
expect the kinetic inductance fraction a to have the same geometrical scaling as g m . 

Fig. 3 shows that the four quantities, F*, a, g m and g g , all scale as s' 1 with 7 = 0.85—0.91. 
The finite thickness of the superconducting film is responsible for the deviations from 7=1. 
This is very strong evidence that the TLS have a surface distribution and are not uniformly 
distributed in the bulk substrate. Our data cannot discriminate between a TLS distribution 
on the metal surface from a TLS distribution on the exposed substrate surface (the gap), 
because the corresponding theoretical predictions (g m and g g ) are very similar and both 
agree with the data. Future measurements of resonators with various center strip to gap 
ratios may allow these two TLS distributions to be separated. 

The absolute values of F*S are also of interest. Assuming a typical value of 5 ~ 10~ 2 for 
the TLS-loaded material HI], the measured value of F*5 = 3 x 10~ 5 for the 3 /zm resonator 
yields a filling factor of F* ~ 0.3%. Numerical calculations show that this is consistent with 
a ~ 2 nm layer of the TLS-loaded material on the metal surface or a ~ 3 nm layer on the 
gap surface, suggesting that native oxides or adsorbed layers may be the TLS host material. 

In summary, the anomalous low-temperature frequency shifts of our superconducting 
CPW resonators are well explained by a mode 
face of the CPW. The excess frequency noise 



in which TLS are distributed on the sur- 

y, 



151 ] also displays a strong geometrical 
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dependence: for a fixed internal power, we find that the noise scales as consistent 
with a surface distribution of TLS fluctuators[22|. It therefore seems very likely that the 
TLS causing the anomalous low-temperature frequency shifts are also responsible for the 
excess frequency noise. The use of optimized geometries or non-oxidizing materials may 
therefore offer a route to more sensitive photon detectors. The TLS should also affect the 
resonator dissipation. At the relatively high power levels used in our experiments, TLS dis- 



sipation is strongly saturated [111. 114 and rather high values of the resonator quality factor 
(Q r > 10 5 ) are routinely obtained^ |. However, at low enough microwave power, the TLS 
response should become unsaturated [231] at which point the quality factor should be limited 
to Q r ~ 1/ FS, or around 3 x 10 4 for the 3 /im resonator described in this paper. This low- 
power regime is of direct relevance for quantum experiments [lo| . in which the microwave 
excitation of the resonator consists of one or a few photons. 
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TABLE 1: Values and ratios 



[urn] 


A/ r (4.2 K) 
[MHz] 
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F*5 
xlO" 5 




9m 


Sg 




1 3(jm 


3m,3^ra 


$g,3/^m 


3 /im 


11.1 


1 


2.98 ±0.12 


1 


1 


1 


5 /xm 


7.41 


0.67 


2.00 ±0.07 


0.67 


0.62 


0.64 


10 fim 


4.15 


0.37 


1.10 ±0.03 


0.37 


0.33 


0.35 


20 urn 


2.28 


0.21 


0.54 ± 0.03 


0.18 


0.17 


0.19 


50 fim 


1.02 


0.092 0.24 ± 0.02 


0.08 


0.075 


0.086 
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FIGURE CAPTIONS 



Figure 1 

(Color online) An illustration of the CPW coupler and resonator. The inset shows a cross- 
sectional view of the CPW. The contour of the metal surface and the contour of the exposed 
surface of the substrate are indicated by the solid line and the dashed line, respectively. 

Figure 2 

(Color online) Fractional frequency shift Af r /f r as a function of temperature. Af r /f r is 
calculated using Af r / f r = [f r (T) — / r (800 mK)]// r (800 mK). The temperature sweep is 
in steps of 50 mK from 100 mK to 600 mK, and in steps of 100 mK above 600 mK. The 
markers represent different resonator geometries, as indicated by the values of the center 
strip width s r in the legend. The dashed lines indicate fits to the TLS theory (Eq. [2]). 

Figure 3 

(Color online) The scaling of the measured values of the kinetic inductance a and TLS filling 
factor F*, as well as the calculated values of the CPW geometrical factors g m and g s , are 
shown, as a function of the resonator center strip width s r . The top panel shows the ratios 
of a (ri), F* (r 2 ), g m (r 3 ) and g g (r 4 ) to their values for the 3 /xm resonators. The bottom 
panel shows these ratios normalized by the kinetic inductance ratio r\. 
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Figure 1 
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Figure 3 
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